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We study quantum interference effects on the transition strength for strongly localized electrons 
hopping on two-dimensional (2D) square and three-dimensional (3D) cubic lattices in the presence 
of a magnetic field. These effects arise from the interference between phase factors associated with 
different electron paths connecting two distinct sites. For electrons confined on a square lattice, 
with and without disorder, we obtain closed-form expressions for the tunneling probability, which 
determines the conductivity, between two arbitrary sites by exactly summing the corresponding 
phase factors of all forward-scattering paths connecting them. By analytically summing paths which 
allow backward excursions in the forward-scattering direction, we find that the interference patterns 
between the dominant winding paths are not drastically different from those between the directed 
paths. An analytic field- dependent expression, valid in any dimension, for the magnetoconductance 
(MC) is derived. A positive MC is clearly observed when turning on the magnetic field. In 2D, 
when the strength of B reaches a certain value, which is inversely proportional to twice the hopping 
length, the MC is increased by a factor of two compared to that at zero field. The periodicity in the 
flux of the MC is found to be equal to the superconducting flux quantum hc/2e. We also investigate 
transport on the much less-studied and experimentally important 3D cubic lattice case, where it is 
shown how the interference patterns and the small-field behavior of the MC vary according to the 
orientation of the applied field B. At very small fields, for two sites diagonally separated a distance 
r, we find that the MC behaves as: rB in quasi-lD systems, r 3 ^ 2 B in 2D with B = (0,0, B), and 
rB [r 3 / 2 B] in 3D with B parallel [perpendicular] to the (1, 1, 1) direction. Furthermore, for a 3D 
sample, the effect on the low-flux MC due to the randomness of the angles between the hopping 
direction and the orientation of B is examined analytically. 

PACS numbers: 72.20.Dp, 72.10.Bg 



I. INTRODUCTION 

Electrons moving on a lattice immersed in a magnetic field have attracted much attention due to their relevance to 
many physical problems. In particular, quantum interference (QI) effects between different electron paths in disordered 
electron systems have been a subject of intense study because they play an important role in quantum transport. For 
instance, the QI of closed loops and their time-reversed paths is central to weak-localization phenomena.^! Indeed, 
during the past decade and half, many fascinating phenomena — including universal conductance fluctuations as well as 
magnetic-field and spin-orbit scattering effects on the conductivity — observed in the weakly localized, metallic regime 
have been understood in terms of the QI between different Feynman diffusive paths in backscattering loops (i.e., paths 
bringing an electron back_to_the starting point). Recently, interest has grown in the effects of a magnetic field on 
strongly localized electronsncl with variable-range hopping (VRH) where striking QI phenomena has been observed in 
mesoscopic and macroscopic insulators or strongly disordered compounds: anomalous magnetoresistance, pronounced 
conductance fluctuations^ Aharonov-Bohm oscillations with periods of hc/e and hc/2e, and the Hall effect. This 
strongly localized regimeBcl is less well-understood than the weak-localization case. 

In the strongly localized regime, the major mechanism for transport is thermally activated hopping between the 
localized sites. In the VRH introduced by Mott,Ej localized electrons, whose wavefunctions decay exponentially with 
a localization length £, hop a distance which is many times larger than £. As a result of the balance between the 
probabilities for hopping and thermal activation, Mott derived that in d dimensions the hopping length changes with 
temperature as £(To/T) 1 /( d+1 \ where To is a characteristic temperature. Therefore, the lower the temperature is, the 
further away the electron tunnels in order to find a localized site of closer energy. 
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According to the "critical path analysis"llil arguments, the conductance of the sample is governed by one critical 
hopping event. During this critical phonon-assisted tunneling process, the electron traverses many other impurities 
since the hopping length is very large at low temperatures. While encountering these intermediate scatterers, the 
electron preserves its phase memory. This elastic multiple-scattering is the origin of the QI effects associated with 
a single hopping event between the initial (i) and final (/) sites. The overall tunneling amplitude Tf between the 
sites i and / is therefore the sum of the contributions from all possible paths connecting them.Elu In other words, the 
tunneling probability of one distant hop is determined by the interference of many electron paths between i and /. 
This leads to Mott's law for the temperature dependence of the conductivity in d dimensions £3 

a(T) ~ \T lf \ 2 exp 

It is worthwhile to note that, in the limit of strong localization, the dominant contribution to Tf comes from 
the shortest paths between i and / (i.e., the "directed path model"). In other words, only interference between 
forward-scattering paths need to be taken into account. This is in contrast with weak localization which results from 
backscattering processes on closed paths. The focus of this paper is on the QI effects on Tf and relevant physical 
quantities due to the presence of an external magnetic field. We will use the model proposed in Ref. 3, which is used 
in most of the recent theoretical work in this area. In this model, the impurities are arranged on a regular square 
(cubic) lattice in 2D (3D) and a nearest-neighbor tight-binding Anderson Hamiltonian is employed. 

In this work we investigate the QI of strongly localized electrons by doing exact summations over all forward- 
scattering paths between two arbitrary sites. We derive compact closed-form expressions for various physical quantities 
(e.g., the transition strength which determines the conductivity) for electrons propagating on a square lattice subject 
to an external magnetic field, with and without random impurities. We also obtain an explicit formula for an 
experimentally important case that has been much less studied theoretically so far: the interference between paths 
on a 3D cubic lattice. 

In the disordered case, by analytically computing the moments for the tunneling probability and employing the 
replica method, we derive analytic results for the magnetoconductance (MC) in terms of sums-over-paths, which are 
applicable in any dimension. Our explicit field-dependent expressions for the MC provide a precise description of the 
MC in terms of the magnetic flux. A positive MC, with a saturation value slightly larger than twice the MC at zero 
field, is observed when turning on the field B. In 2D, the saturated value of the magnetic field -B sa (i.e., the first 
field that makes the MC become twice the value at zero field) is inversely proportional to twice the hopping length: 
the larger the system is, the smaller B sa will be. In other words, as soon as the system, with hopping distance r, is 
penetrated by a total flux of (r/8)(hc/e), the MC reaches the saturation value. The period of oscillation of the MC 
is found to be equal to hc/2e, which is the superconducting flux quantum. 

Furthermore, at very small fields, for two sites diagonally separated a distance r, the MC scales as follows: (i) r B 
for quasi-lD ladder-type geometries with B = (0,0, B), (ii) r 3 / 2 B in 2D with B = (0,0,5), (iii) rB in 3D with B 
parallel to the (1, 1, 1) direction, and (iv) r 3 / 2 B in 3D with B perpendicular to the (1, 1, 1) direction. 

The general expressions presented here: (i) contain, as particular cases, several QI resultsou derived during the 
past decade (often by using either numerical or approximate methods), (ii) include QI to arbitrary points (m,n) on 
a square lattice, instead of only diagonal sites (m,m), (iii) focus on 2D and 3D lattices, and (iv) can be extended to 
also include backward excursions (e.g., side windings) in the directed paths. 

Exact results in this class of directed-paths problems are valuable and can be useful when studying other systems, 
for instance: (1) directed polymers in a disordered substrateisee, e.g., Refs. 12 and 13), (2) interfaces in 2D (see, for 
instance, Ref. 14), (3) light propagation in random media,li§E3 and (4) charged bosons in 1D.0 

To study the magnetic field effects on the tunneling probability of strongly localized electrons, we start from the 
tight-binding Hamiltonian 

H = WJ2 4* + V c\c, exp(z^) , (1) 

« (*J> 

where (ij) refers to the nearest-neighbor sites and the phase Ay = 2tt J j A-dl is 2ir times the line integral of the 
vector potential along the bond from i to j in unitji..of the normal flux quantum hc/e. In the strongly localized regime, 
V/W <C 1, the electron energy can be set to zero.tau Consider two states, localized at sites i and / which are r bonds 
apart. By using a locator expansion, the transition amplitude (Green's function) Tf between these two states can be 
expressed ascltl 

°° f V \ r+21 

T «=Y. w (w) s(r+2l) > (2) 
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where 



S (r+2l) = J- e** r , (3) 

All (r+2/)-step paths T 
connecting sites i and / 

and $r is the sum over phases of the bonds on the path T of r + 2l steps connecting sites i and /. In general, T contains 
closed loops. In the strongly localized regime (i.e., V/W <C 1), the dominant contribution to Tif is W(V/W) r S^ r \ 
where 

Sir) = e l * r , (4) 

All directed paths T 
of r steps on a lattice 

In other words, only the shortest-length paths (with no backward excursions) connecting them arc taken into account, 



lr-model provides an excellent approximation to j£ 
£j It is important to stress that the conductivity^ 



namely, the directed-path model of Refs. 2-8. This directed-pa 
since (V/W) 2 is quite small in the extremely localized regime, 
between i and / is proportional to |Tiy| 2 . 

Quantum interference, contained in S^ r+2l \ arises because the phase factors of different paths connecting the 
initial and final sites interfere with each other. We will first focus on the computation of S^, which is the essential 
QI quantity for electrons deep in the localized regime. In 2D, we also analytically compute S , ( r + 2 ) which becomes 
important when electrons are not so strongly localized. 

This paper is organized as follows. In Sec. H,r5jv e study QI on a square lattice under a uniform potential, which is 
related to the decay of gap states into the bulk.El Here we derive an elegant, general, and very compact closed-form 
expression for S^ r \ Intriguing properties associated with the behavior of on diagonal sites are discussed in detail. 
It will be shown later (in section III) that the effect of a magnetic field on the MC is governed by the behavior of S^'K 
As a step towards the understanding of interference between non-directed paths, we also go beyond the directed-path 
model by exactly computing analytic results for S^ r+2 \ 

In Sec. Ill, we investigate the tunneling in a random impurity potential, which is relevant to the conductivity of, 
for example, lightly doped semiconductors and strongly disordered compounds.cl Closed-form results for the tunneling 
probability, which determines the conductivity, are obtained. We then analytically compute the moments for the 
tunneling probability. From them, we derive analytic field-dependent expressions, valid in any dimension, for the 
MC. The full behavior of the MC as a function of the magnetic flux — including the scaling in the low-field limit and 
the occurrence of saturation — is discussed in detail. The close relationship between the QI quantity and the 
corresponding MC is illustrated. Comparison of our results with experimental observation and other theoretical work 
is also made. 

In Sec. IV, we examine the QI on a 3D cubic lattice and provide a general formula for S^ r \ We show how the 
interference patterns and the small-field behavior of the MC vary according to the orientation of the applied field. 
Furthermore, we investigate the effect on the low-flux MC due to the randomness the angles between the directions 
of the critical hops and the orientation of the applied field. 

In Sec. V, we conclude by addressing several relevant issues and summarize our results. 



II. QUANTUM INTERFERENCE ON A TWO-DIMENSIONAL SQUARE LATTICE 

A. Exact summation of forward-scattering paths: S^ r ' 

Let (m, n) denote the site coordinates. Without loss of generality, we choose (0, 0) to be the initial site and focus on 
m,n > 0. For forward-scattering paths of r steps, which exclude backward excursions (i.e., only moving upward and 
to the right is allowed), ending sites (m, n) satisfy m + n = r. Let S m n (—S^) be the sum over all directed paths of 
r steps on which an electron can hop from the origin to the site (m, n), each one weighted by its corresponding phase 
factor. Employing the symmetric gauge 

A = yH/>a), 

and denoting the flux through an elementary plaquette (i.e., with an area corresponding to the square of the average 
distance, which is typically equal to or larger than the localization length £, between two impurities) by 4>/2t:, it is 
straightforward to construct the recursion relation: 
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c _ -in<t>/2 o i im4>/2 c (k\ 

^m,n — c ^m— l,n ' c ^m^n—l • \° J 

This equation states that the site (m, n) can be reached by taking the rth step from neighboring sites to the left or 
below. The factors in front of the S"s, namely, exp(— incfi/2) and exp(irn</>/2), account for the presence of the magnetic 
field. Enumerating the recursion relations for Sk n , n (k n = fn — 1, m — 2, . . . , 0) successively and using So, n = 1 for any 
n, we obtain the following relation 

m 

S m ,n = e-^ m ~ k ^' 2 S kn , n -l . (6) 

fc„=0 

Here SVn.n is expressed as a sum of the S"s one row below (i.e., on the line y = n— 1). The physical meaning of Eq. (6) 
is clear: the site (m, n) can be reached by moving one step upward from sites (k n , n — 1) with < k n < m, acquiring 
the phase ik n <fi/2; then traversing m — k n steps from (k n ,n) to (m, n), each step with a phase —inip/2. By applying 
Eq. (6) recursively (and utilizing S m ,o = 1 for any m), 5 m ,n for m, n > 1 can be written as 



iW = e --« ^/ 2 i mi „(0), (7) 



where 



L m .n(4>)=Y, E ' ' ' E e l ( fel + - + fe "- 1+fc "^ (8) 
fc n =0fe n _i=0 fei=0 
n fcj + i 

=n e 

j=l /c 3 =0 

with k n -\-\ = m. If we use the Landau gauge instead, the expression for the sum-over-paths will read L m/n ] 
namely, S m>n employs the Symmetric gauge, while L m>n uses the Landau gauge. Notice that each term in the 
summand corresponds to the overall phase factor associated with a directed path. In the absence of the magnetic flux 
(0-0), 

Sm,n(0) = E E •••E 1 
k n — k n —i= ki=0 

= c™ +n = = N, (9) 

m! n\ 

which is just the total number of r-step paths between (0, 0) and (m, n). 

After some calculations we obtain one of our main results, a very compact and elegant closed- form expression for 

S m ,n{4>)'- 

*"<*>- (10) 



where 



m , 

F m {4>)= II sin o0- (11) 

fc=l 1 

Notice that the symmetry S m . n — S n , m [apparent in Eq. (10)] is due to the square lattice geometry. Moreover, we 
also obtain 

YYk-i (i _ eifc0 ) 

= [nr =1 (i-e**)][nLi(i-e ifc *)] ■ (12) 

Previous work on QI in the VRH regime obtained particular cases, mostly numerical, of sums to diagonal points 
S m ,m, while the general result Eq. (10) is valid for arbitrary (e.g., non-diagonal) sites. 

To illustrate the quantum interference originating from sums over phase factors associated with directed paths, we 
show the five possible ending sites (m, n) for r — 4 and their corresponding SVn.n in Fig. 1(a). In addition, the six 
different paths connecting (0, 0) and (2, 2) and their separate phase factor contributions to are shown in Fig. 1(b). 
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B. Low- flux limit 



In the very-low-flux limit <gc 1, the logarithm of 5 m; „, calculated exactly to order cj) 2 , is 

lniS m)n (^i) = InN — —mn (m + n + 1) <j) 2 . (13) 

Thus we obtain the familiarH harmonic shrinkage of the wave function with explicit expressions for all the prefactors. 
This result can be interpreted as follows. The effective "cigar-shape" area exposed to the field has an effective length 

left' 

l e S ~ y/mn 

(i.e., the square root of the area enclosed by the paths) and an effective width w e g: 



WeE ~ \/m + n 

(i.e., the square root of the length of the paths). For the special case m = n = r/2, 

lnS m , m (<t>) = In Jl - ^r 2 (r + 1) 4> 2 , (14) 

which is consistent with, and generalizes, the results in Ref. 4 since it gives the exact prefactor. Thus, the effective 
length is ~ r, while the effective width is ~ yfr. Furthermore, for a ladder-type quasi-lD system, (e.g., m = r — 1 
and n = 1), we have 

lnS r _i,i($ =lnr-^(r 2 -l)^ 2 . (15) 

In this case, the effective length and width are both ~ y/r. This result remains valid for small values of n (narrow 
stripes or multiladders) . The fourth-order contribution to In 5 mjrl (</>) can also be computed exactly as 

[ 10 (6m 4 + 15m 3 n + 20m 2 n 2 + 15mn 3 + 6n 4 ) 
+ 114 (to 3 + 2m 2 n + 2mn 2 + n 3 ) 
+ 29 (2m 2 + 3mn + 2n 2 ) + 9 (m + n) + 5 ] 4 . 



C. Quantum interference on diagonal sites 

Among the S"s for an even number of steps, those located along the diagonal corners contain the richest interference 
effects since the number of paths ending at (to, to) and the areas they enclose are both the largest. We therefore 
examine more closely the behavior of the quantities 

i 2m {4>) = s m , m {4>) = n sm ~r . (i6) 

For irrational flux 0, it can be proved that — 1 < li m < 1 for any m. A particular case (asymptotic behavior) of our 
very compact general expression Eq. (16) for l2 m is investigated in detail by Fishman, Shapir, and Wang in Ref. 6. 
For <j) = 2tt s/t (s and t are positive integers with 1 < s < t and s being prime to t), we obtain (n > 0) for to <t 



_ (-l)""(2n)! T 
J-K m.4-nt\ — : — : J; 

and 

^,,=(-1)*"'^. ( I S ! 

Furthermore, for those to satisfying 



2(m+nt) — i — j Hmi \± I ) 



(nl) 2 ' 
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^<m<t-l, (19) 

h(m+nt) = • (20) 

In other words, the zeros of l2m{4>) are given by <fi = 2tt s/t for 

at±l< t <-^-, (2i) 

n + 1 ~ _ In + 1 ' v 7 

with 

777 — 1 

< n < HL- , (22) 

and the s's are prime to each allowed t. From a physical viewpoint, these flux values produce the complete cancelation 
of all phase factors (i.e., fully destructive interference) and result in the vanishing of the tunneling probability (and 
conductivity). Indeed, we will see in section III that if we also consider the effects of the on-site impurity scattering, 
these flux values lead to the largest (i.e., saturated) value for the positive MC. 

In Fig. 2, we show the zeros for I2, I4, Iq, . . . , I40 obtained by using Eqs. (21) and (22). Note that the smallest value 
of 4>/2ir satisfying l2m{4>) = is always l/2m and the number of zeros rapidly increases when m becomes larger. 

The 12m can be expressed as sums of trigonometric cosines. For instance, the first few ones are (with = <t>/2): 

h = 2 cos e, 

I 4 = 2 + 2cos20 + 2cos40, 

I e = 6cos0 + 6cos30 + 4cos50 + 2cos70 + 2cos90, 

I s = 8 + 14 cos 29 + 14 cos 40 + 10 cos 60 + 10 cos 86 + 6 cos 100 + 4 cos 126 + 2 cos U6 + 2 cos 166, 
I 10 = 40 cos 6 + 38 cos 36 + 36 cos 56 + 32 cos 76 + 28 cos 96 + 22 cos 116+ 18 cos 136 

+ 14 cos 156* + 10 cos 176* + 6 cos 196* + 4 cos 216* + 2 cos 236* + 2 cos 256, 
J12 = 58 + 110 cos 20 + 110 cos 46 + 102 cos 66 + 96 cos 86* + 84 cos 106* + 78 cos 120 + 64 cos 140 

+ 56 cos 160 + 44 cos 180 + 36 cos 200 + 26 cos 220 + 22 cos 240 + 14 cos 260 

+ 10 cos 280 + 6 cos 300 + 4 cos 320 + 2 cos 340 + 2 cos 360, 
I 14 = 338 cos + 332 cos 30 + 324 cos 50 + 310 cos 70 + 292 cos 90 + 272 cos 110 + 250 cos 130 + 224 cos 150 

+ 200 cos 170+ 174 cos 190 + 150 cos 210 + 126 cos 230 + 106 cos 250 + 84 cos 270 + 68 cos 290 + 52 cos 310 

+ 40 cos 330 + 30 cos 350 + 22 cos 370+ 14 cos 390 + 10 cos 410 + 6 cos 430 + 4 cos 450 + 2 cos 470 + 2 cos 490. 

Notice that hm depends only on the even (odd) harmonics of when m is even (odd). I2m{4>) obeys the following 
properties: (i) 2tt (4tt) periodicity in <f> for even (odd) to, namely, 

/ 4 ™(0 + 27r) =I 4n (<P), 

hn+l{4> + 471") = hn+2{4>) ■ 

In other words, the period of /2m corresponds to hc/e when to is even, and 2hc/e when to is odd. (ii) With to even 

hm(2n -<p) = h m {4>) 

for < 4> < 7r. Also 

to! 



[(tow 

(iii) With to odd 

/ 2m (27T±0) = -I 2m (<j>) 

for < <p < 2ir. 

From the properties described above, we can draw a general picture of the behavior of hm- Let $ = 4>/2n. 
I2m(& = 0) = (2to)!/(to!) 2 , which is a enormous number for large to. As the magnetic field is turned on, /2m rapidly 
drops to its first zero at <f> — l/2m. I 2m then shows distinct behaviors depending on m. 
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For even m, hm exhibits many small-magnitude fluctuations around zero for =^ < $ < ™ ■ li m then monotoni- 
cally climbs from to a large positive value, hm{^) = to!/[(to/2)!] 2 , for < $ < i. It is evident that /2m(7r) is 
still very small compared to /2m (0). Within the period < $ < 1 (i.e., < 4> < 2tt), I 2m has mirror symmetry with 
respect to $ = 1/2 (i.e., cf> — it). 

For odd to, li m exhibits many small-magnitude fluctuations around zero for -L- < $ < 2 " m 1 . In addition, 7 2 m 
always equals at $ = 1/2 (i.e., <j> = tt). hm then monotonically drops from to -(2m)!/(m!) 2 for 2 ™~ 1 < $ < 1. 
For < 4> < 1, /2m has inversion symmetry with respect to $ = 1/2. Within the period < $ < 2 (i.e., < <fi < 47r), 
J 2m has mirror symmetry around $ = 1 (i.e., <j) = 27r). Recall that for any to, — 1 < /2m < 1, for irrational values of 

In Fig. 3, we plot I 2 through lis, ho, hs an d ^40- These figures show very interesting interference patterns 
of J 2m and clearly reflect the general description given above. It is worthwhile to keep in mind that the properties 
embedded in S m ,n described above play a central role in determining the behavior of the MC obtained in section III. 



D. Exact summation of the dominant winding paths: 5'' r+2 ' 

Up to now, we have focused on the computation of and presented a detailed investigation of their properties. 
When the electrons are less strongly localized, the next higher-order contribution to Tif (i.e., W(V/W) r+2 S^' +2 \ 
which is the dominant term including backward excursions) becomes important. Therefore, quantum interference 
effects between phase factors of paths with backward recursions (i.e., moving downward and to the left is also included) 
need to be taken into consideration. Notice that paths in S^ r+2 \ though include backscattering processes, do not 
involve closed loops enclosing flux (e.g., elementary square plaquettes). 

In this section we present the computation of the second-order contribution, namely S^- r+2 \ to the transition 
amplitude T^-. Let P m ,„ (= S^- r+2 ^) denote the sums over paths of to + n + 2 steps starting from (0,0) and ending 
at (m,n). We assume that the electrons are confined on a square lattice with non-negative x and y coordinates. We 
can divide the contribution to P m ^ n into five parts. 

First, hopping directly to site (p, 0), with 1 < p < to, electrons take one step back to (p— 1, 0), then hop to— p+ l + n 
steps to (to, n). Second, hopping directly to site (0, q), with 1 < q < n, electrons take one step back to (0, q — 1), then 
hop to + n — q + 1 steps to (to, n). Third, directly hopping to site (p, n + 1), with < p < to, electrons move one-step 
downward to (p, n) gaining a phase factor exp(—ip<f)/2), then hop to — p steps to (m,n). Fourth, directly hopping 
to site (to + l,q), with < q < n, electrons move one-step downward to (m,q) gaining a phase factor cxp(iq4> / 2) , 
then hop n — q steps to (to, n). Fifth, directly hopping to (p, q) with 1 < p < m and 1 < q < n, electrons take one 
step back to either (p — l,q) or {p, q — 1), accompanied by the phase factor exp{iq4>/2) or exp(— ip<j>/2), then hop 
TO + n — p — q + 1 steps to the ending site (to, n). Therefore P m ,n can be written as 

m n m n 

p—1 q—1 p—0 q—0 



m n 



~l~ ^ ' ^""^ Sp.q |f Bp—l q^ rn n +e Bp q—l^ m „) . (23) 

P=l 9=1 

where B VA -, m ^ n is the sum over phase factors of all directed paths (i.e., containing to + n — p — q steps) starting from 
(p, q) and ending at (to, n). After some calculation we obtain 

R \ A-{m-p)q+{n- q)p](t> \ 

^p.q^-m.n — LA[J s t ^ ? o m —p.n — q- V ^7 

By substituting Eq. (24) into Eq. (23), we derive 

m n m n 

V in( V -l)4> v _ i m (q-l)4> _ imrL<P v i(^-l)p</> imnj, v _ i(m-l)q<p 

Pra.n = e 2 S m -p+l,n + }^ G * Sm,n-q+l + € 2 ^ e 2 S P,n+l + e 2 e 2 S m+ i. q 

p—1 q—1 p—0 q—0 

m n 

+ E E S ™ {e-^ m -^-^v-m/2 Sm _ p+1 , n _ q + e -iM*-U-(»-Ui#/2 Sm -p, n - q+1 } . (25) 

P=l 9=1 

In the special case to = n, 



7 



„ \m(m — j\ 

2 E C0S 2 
i=i L 

m — 1 J — 1 

~l~ 2 ^ ^ ^ ^ $m—j,m—k 
j=l k=0 



COS 



+ 2 S*\ cos 



(TO 2 — TO J + j)0 



+ 2 cos 



2 I >->m-j,m-3 



(mj — mk — j) 



~>j,k+l + COS 



(mj — mk + k)4> 



Sk 



(26) 



The explicit expressions for the first few P m%m are (with = <j)/2): 



4 cos 130 + 2 cos 15(9, 
- 130 cos 126» 



14 cos + 2 cos 36», 

26 + 32 cos 26* + 26 cos 40 + 4 cos 60 + 2 cos 80, 

130 cos + 124 cos 30 + 88 cos 50 + 52 cos 70 + 40 cos 90 + 8 cos 110 4 
224 + 410 cos 20 + 396 cos 40 + 308 cos 60 + 282 cos 80+ 188 cos 100 
+ 76 cos 140 + 58 cos 160 + 14 cos 180 + 8 cos 200 + 4 cos 220 + 2 cos 240, 
1446 cos + 1386 cos 30 + 1308 cos 50 + 1176 cos 70 + 1032 cos 90 + 842 cos 110 + 690 cos 130 
+ 542 cos 150 + 398 cos 170 + 264 cos 190 + 180 cos 210 + 108 cos 230 + 80 cos 250 
+ 24 cos 270 + 14 cos 290 + 8 cos 310 + 4 cos 330 + 2 cos 350, 

2518 + 4868 cos 20 + 4808 cos 40 + 4514 cos 60 + 4238 cos 80 + 3788 cos 100 + 3466 cos 120 
+ 2938 cos 140 + 2554 cos 160 + 2074 cos 180 + 1702 cos 200 + 1298 cos 220 + 1056 cos 240 
+ 736 cos 260 + 536 cos 280 + 356 cos 300 + 244 cos 320+ 148 cos 340 + 1 10 cos 360 
+ 38 cos 380 + 24 cos 400+ 14 cos 420 + 8 cos 440 + 4 cos 460 + 2 cos 480. 

These -P m , m 's are plotted in Fig. 4. Note that P m ,m depends only on the even (odd) harmonics of and has a period 
27r (47r) for even (odd) to. The expressions for P m . m are obviously more complicated than the corresponding 12m- 
However, by comparing Figs. 3 and 4, we find that the general features in the interference behaviors are surprisingly 
similar. We thus infer that the relevant physical quantities are not significantly changed by the addition of interference 
between the dominant winding paths. 



Phi 

-P3,3 
Pi,i 

P5,5 



6,6 



III. EFFECTS OF DISORDER 



A. Average of the tunneling probability 

To incorporate the effects of random impurities, we now replace the on-site energy part in Eq. (1) (first term in H) 
by J2i e i c l c i> where the q's are now independent random variables. The Hamiltonian now takes the form 

H = ^2 e i C i C i + V Y1 C * C J e *P( iA v) ■ 

We have studied two commonly used models: (i) a can take two values: +W and — W with equal probability; and 
(ii) a is randomly chosen from a uniform distribution of width W and zero mean. We found that both models yield 
the same results for the MC. 

We now start with the general case of the first model, namely, can have two values: +W with probability fi and 
— W with probability v, where [i + v = 1. Due to disorder, the transition amplitude becomes 

T if = w(^j J m , n , 

with 

(27) 

where T runs over all directed paths of r steps connecting sites (0,0) and (m, n), and j over sites on each path. For 
all directed paths ending at (m,n), electrons traverse r = to + n sites (the initial site (0,0) is excluded). Each site 




8 



visited now contributes an additional multiplicative factor of cither +1 or —1 to the phase factor. Therefore, for a 
given path T, the probability for obtaining ±e l * r is 

= ^ + uY±^-vY 

± 2 

It is then clear that 

( J m , n (4>) ) = (P+- P-) Sm,n(<f>) = (P- "Y S m , n (4>), (28) 

where (• • •) denotes averaging over all impurities. 

By exploiting Eqs. (7) and (8), we derive the following general expressions valid for any fj, and v, 

( Jl.M ) = (1 - V) S m , n (2<j>) + V S 2 m .M), (29) 

where 

JV-l 

r = pN + pN + Y^ pN-k pk {C N _ 4 c N_- 2) 
fe=l 

= 1 -4P+P_ = ( M -z/) 2r . (30) 
Also, the disorder average of the tunneling probability (i.e., the transmission rate) \ J\ 2 = J J* is 

( \JmA4>)\ 2 ) = (1 - V) N + V SlJcj,). (31) 
The physical origin of Eqs. (29) and (31) becomes clearer by rewriting them as 

( ) = S rnA 2 <P) + T> Sm,n(<P) [S m ,n{<l>) - C m ,n(4>)} (32) 

and 

(\J m ,n(<l>)\ 2 ) = N + V[S 2 m>n (<j>) - N], (33) 

where 

ro '" w "(nL 1 cos|<A)(nL 1 cos|<A)' 

and we have used 

S m , n (2<f>) = S m . n ((t>)C m . n (4>). 
The first terms in Eqs. (32) and (33) account for contributions from pairs of identical paths: 

]T (±e** r ) (±e** r ) = £e 2i * r = S m ,n(W) 
r r 

in Um,J> and 

]T(±e 4 * r ) (±e- l * r ) = 1 = N 
r r 

in ( |J m ,ri| 2 )- The second terms in Eqs. (32) and (33) account for contributions from pairs of distinct paths. Note 
that S m ^ n (0) = N and C TOi „(0) = 1, when = 0. We then have in the absence of magnetic flux 

( Jm,n(°) > - ( \Jm, n m 2 > = N + V N (N - 1) . (34) 

Furthermore, in the special case = v = 1/2, since V = we then obtain 

( J m,n{4>)) = S m ,n(2<j>), 

(\J m A4>)\ 2 )=N. 
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B. Higher-order moments and general expressions for the first few leading terms 



For fi = v = 1/2 (the most studied case so far), we can obtain analytical expressions for the moments ( J^.nW) ) 
and ( \ J m ,n{4 l )\ 2p ) f° r an Y value of p. Only a few of these will be presented here. From now on, </(</>) stands for 
Jm,n(4>) an d S((j>) stands for 5 TOj „(0). The derivation of these moments is given in appendix A. 

( J 4 (0) ) = 35 2 (20) - 2 5(40), (35.1) 

( J 6 (0) ) = 15 5 3 (20) - 30 5(20) 5(40) + 16 5(60), (35.2) 

( J 8 (0) ) = 105 5 4 (20) - 420 5 2 (20) 5(40) + 448 5(20) 5(60) + 140 5 2 (40) - 272 S(8(f>), (35.3) 
( J 1O (0) ) = 945 5 5 (20) - 6300 5 3 (20) 5(40) + 10080 5 2 (20) 5(60) + 6300 5(20) 5 2 (40) 

- 12240 5(20) 5(80) - 6720 5(40) 5(60) + 79365(100), (35.4) 
( J 12 (0) ) = 10395 5 6 (20) - 103950 5 4 (20) 5(40) + 221760 5 3 (20) 5(60) + 207900 5 2 (20) 5 2 (40) 

- 403920 5 2 (20) 5(80) - 443520 5(20) 5(40) 5(60) - 46200 5 3 (40) 

+ 523776 5(20) 5(100) + 269280 5(40) 5(80) + 118272 5 2 (60) - 353792 5(120), (35.5) 
( J 14 (0)) = 135135 5 7 (20) - 1891890 5 5 (20) 5(40) + 5045040 5 4 (20) 5(60) + 6306300 5 3 (20)5 2 (40) 

- 12252240 5 3 (20) 5(80) - 20180160 5 2 (20) 5(40) 5(60) - 4204200 5(20) 5 3 (40) 
+ 23831808 5 2 (20) 5(100) + 24504480 5(20) 5(40) 5(80) + 10762752 5(20) 5 2 (60) 
+ 6726720 5 2 (40) 5(60) - 32195072 5(20) 5(120) - 15887872 5(40) 5(100) 

- 13069056 5(60) 5(80) + 22368256 5(140), (35.6) 

and 

(|J(0)| 4 )=2TV(TV-l) + 5 2 (20), (36.1) 

( | J(0)| 6 ) = 2 TV (37V 2 - 97V + 8) + 3 (37V - 4) 5 2 (20), (36.2) 
( | J(0)| 8 ) = 8iV (37V 3 - 187V 2 + 417V - 34) + 8 (97V 2 - 337V + 32) 5 2 (20) 

+ 9 5 4 (20) - 12 5 2 (20) 5(40) + 4 5 2 (40), (36.3) 
( | J(0)| 10 ) = 87V (157V 4 - 1507V 3 + 6257V 2 - 12507V + 992) + 40 (157V 3 - 1057V 2 + 2607V - 216) 5 2 (20) 

+ 75 (37V - 8) 5 4 (20) - 20 (157V - 44) 5 2 (20) 5(40) + 20 (57V - 16) 5 2 (40), (36.4) 
( | J(0)| 12 ) = 167V (457V 5 - 6757V 4 + 44257V 3 - 155257V 2 + 287067V - 2212) 

+ 120 (457V 4 - 5107V 3 + 22957V 2 - 47027V + 3552) 5 2 (20) + 30 (1357V 2 - 8547V + 1440) 5 4 (20) 
+ 225 5 6 (20) - 120 (457V 2 - 3097V + 556) 5 2 (20) 5(40) - 900 5 4 (20) 5(40) 
+ 24 (757V 2 - 5557V + 1064) 5 2 (40) + 900 5 2 (20) 5 2 (40) + 480 5 3 (20) 5(60) 

- 960 5(20) 5(40) 5(60) + 256 5 2 (60), (36.5) 
( | J(0)| 14 ) = 16 TV (3157V 6 - 66157V 5 + 624757V 4 - 3344257V 3 + 10573227V 2 - 18541607V + 1398016) 

+ 56 (9457V 5 - 15757V 4 + 1107757V 3 - 4017307V 2 + 7325367V - 518464) 5 2 (20) 
+ 1470 (457V 3 - 4957V 2 + 19207V - 2576) 5 4 (20) + 11025 (TV - 4) 5 6 (20) 

- 840 (105TV 3 - 12397V 2 + 50967V - 7192) 5 2 (20) 5(40) - 2940 (157V - 64) 5 4 (20) 5(40) 
+ 2940 (157V - 68) 5 2 (20) 5 2 (40) + 3360 (77V - 31) 5 3 (20) 5(60) 

+ 56 (5257V 3 - 66157V 2 + 287847V - 42688) 5 2 (40) - 448 (1057V - 493) 5(20) 5(40) 5(60) 

+ 1792 (77V -34) 5 2 (60). (36.6) 

These moments satisfy the consistency check ( J 2p (0) ) = ( | J(0)| 2p ) and odd moments vanish by symmetry. 

The moments provide an analytical view of the structure of the QI in the tunneling process. Since \ J\ 2 = J J*, 
each \ J\ 2 represents TV forward paths to (m,n), each one with its corresponding reversed path back to the origin. 
Also, ( |J(0)| 2p ), averages over the contributions of TV P such pairs of paths. In general, ( |J(0)| 2p ) consist of terms 
involving N k (k = 1, • • • ,p). The above explicit expressions for the moments will allow us to deduce general formulae 
for the first few leading (i.e., dominant) terms in the moments. 

We first focus on the leading terms (oc TV P ), since they provide the most significant contribution to the moments 
when TV is large. Recall that 5(0) = TV, therefore we need to consider all terms involving 5 2fc (20) N p ~ 2k in ( | J(0)| 2p ). 
We derive (see appendix B for more details) 
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(|J(0)n = (2p-l)UN*>, 



2 A 



, k=0 



S(24>) 



N 



2k' 



Furthermore, by considering all the second leading terms (oc N p 1 ), we obtain 
1 

3' 



\j(<p)\ 2p ) 



6 



E 

,fc=0 



(2fc + l)! 



(2 fc fc!) 2 



(3p + 2fc - 3) Cf fe+1 + 2fcC* 



5(40) 



2fc+2 



TV 



5(20) 



n 2fc 



N 



(37) 
(38) 

(39) 
(40) 



Also, when 5(2p0) — with p > 1 (e.g., at = 7r/m, S m . m (2p(f>) — 0), the third leading terms (oc N p 2 ) in the 
moments are 



( \A0)\ 2p ) = (p - 1) (p - 2) (5p + 1) [(2p - 1)!!] A^ 2 , 
( 1^) I 2p ) = -^P (P ~ 1) (p - 2) (9p + 5) (p!) 7V^ 2 . 



(41) 
(42) 



The above general expressions for the moments are of value since they enable us to analytically obtain the dominant 
contributions to the quantity we are interested in: the magnetoconductance. 



C. Analytical results for the magnetoconductance 

We now use the replica method: 



(ln|J(0)| 2 ) = lim 



|J(0)| 2p )-l 



p^o p 

to compute the log-averaged MC with respect to the zero-field log-averaged MC (denoted by Lmc), defined as 

L MC ee (In | J(0)| 2 )- (In | J(0)| 2 ) 
;|J(0)| 2 P)-(|J(O)| 2 ^) 



(43) 



lim ■ 



P 



(44) 



Taking into account only the first leading terms in the moments, shown in Eqs. (37) and (38), we derive the Lmc as 



, o ^ (2fc-l)! 

L MC = m2-^72W- 

fe=i v ; 



5(20) 
N 



2k 



where we have 



fSJ.], : . 



Exploiting the following identity^ for 



k=l 



< x < 1 



(45) 



(46) 



cosh 1 — = In ^7-7- 



(2k-iy. x2k 



r a j-^ (2 k k\) 



which reduces to Eq. (46) for x = 1, we thus obtain a very concise exact expression for the Lmc as 



Lmc — 



cosh 
In 2 



-1 AT 



\S(24>)\ |S(20)| 



ln ^m when 5(20)^0 
when 5(20) = 



In 



(47) 
(48) 
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The typical MC of a sample 

G(0)=cxp((ln|J(0)| 2 )) 
is then given by, normalized by the zero-field MC G(0), 



G{4>) 
G(0) 



exp(L M c) = 1 + 



N 



(49) 



Eq. (49) is one of our main results. It provides a concise closed-form expression for the MC, as an explicit function 
of the magnetic flux. From Eq. (49) it becomes evident that a magnetic field leads to an increase in the positive MC: 
G(0)/G(O) increases from 1 to a saturated value 2 (since S(2<j>) decreases from TV to 0) when the flux is turned on 
and increased. G{4>) = 2G(0) at the field 4> that satisfies S(2<f>) — 0. Furthermore, it is clear that the MC varies 
periodically with the magnetic field and the periodicity in the flux is equal to the superconducting flux quantum hc/2e. 

It is important to point out that Eqs. (48) and (49) are valid in any dimension as long as we use the corresponding 
D-dimensional sum . 

It is illuminating to draw attention to the close relationship between the behaviors of I 2m (2(p) = S m , m {24>) and 
the corresponding G(<p). When = 0, (I2m(0)/N) 2 = 1, which is the largest value of (I 2m (2(j))/N) 2 as a function of 
(f>, and the MC is equal to the smallest value G(0). When the magnetic field is increased from zero, (I 2m (2(f)) /N) 2 
quickly approaches (more rapidly as m becomes larger) its smallest value, which is zero, at <f>/2i: = l/4m. At the 
same time, the MC rapidly increases to the largest value 2G(0). 

The physical implication of this is clear: fully constructive interference in the case without disorder leads to the 
smallest hopping conduction in the presence of disorder. While fully destructive interference in the case without 
disorder yields the largest hopping conduction in the presence of disorder. Moreover, when m (the system size is 
to x m) is large, G(0)/G(O) remains in the close vicinity of 2 for <P/2t: > l/4m in spite of the strong very-small- 
magnitude fluctuations of I 2m (2(f>)/N around zero. 

The saturated value of the magnetic field B sa (i.e., the first field that makes G(<f>) = 2G(0)) is inversely proportional 
to twice the hopping length: the larger the system is, the smaller B sa will be. In other words, as soon as the system, 
with hopping distance r — 2m, is penetrated by a total flux of (l/2r) x (r/2) 2 = r/8 (in units of the flux quantum 
hc/e), the MC reaches the saturation value 2G(0). 

Defining the relative MC, AG(<f>), as 



AG(0) 



- G(0) 
G(0) 



and utilizing Eq. (49), we show AG((f>) versus <j> for several different hopping lengths in Fig. 5. The behavior of AG(<f>) 
described above can be clearly observed in these figures. 

Now let us examine the behavior of the MC in the low-flux limit. From Eqs. (14) and (15), it follows then that, for 
very small fields, in 2D 



AG(4>) 



^3 r3/2 ( 



(50) 



and in ladder-type quasi-lD structures 



(51) 



In Fig. 6, we plot AG(4>) computed directly from Eq. (49), for various small values of <f>, versus r 3 / 2 0, with r = 
2,4,..., 1000, in (a) and versus rcj), with r — 2,3, ... , 500, in (b), respectively for 2D and quasi-lD systems. It is seen 
that, both in (a) and (b), all the data nicely collapses into a straight line, which verifies the scaling of the low- flux 
MC in Eqs. (50) and (51). 

If we consider the second leading terms in the moments, namely Eqs. (39) and (40), the second-order contribution 
to the Lmc is 



Lmc 



677 j 1 _ Sfc=l 



(2fc+l)! / 2fc-3 k S(44,) \ ( S(24>) \ 

(2 k fc!) 2 ^2fc+l fe+1 N J \ - ' 



2 k 



when S(2<P)/N = ±1 
when S(2<P)/N^±1 



l-A 
A a (l+A) 



N 



k+l N J \ N J 

when A = 
) (1 + 2A) - A 2 (2 + 3A)j } when A ^ ' 



(52) 
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where 



A 



and we have used 



E 



(2fc)! 3 p k _i- VT= 



2 



^ (2 fc fc!) 2 VT 



x 2 



2k 



'( 2 fcfc!)2 (1- ^2)3/2 

(2fc)! - 2 



z 2fe 



E ( fc + i)( 2 fcfc!)2- (i-VT^- . 

The principal features in the behavior of the MC are not significantly modified by the addition of the contribution 
from Lmc in Eq. (52): while the magnitude of AG(cf>) is slightly increased for <fi ^ 0, the period of the MC remains 
unchanged. 

For a 2D system and in the low-flux limit, we derive from Eq. (52) for diagonal sites (r/2,r/2): 

AG(0) = (V3/2AN) r 3/2 <j>. 

Comparing this result with Eq. (50), we see that the dependence of the small- field MC on the hopping length and the 
field is the same except for different prefactors. Summing up both contributions, we have for small <j> 

AG(<t>) = ^ ( 1 + 7^7 V 3/2 *• (53) 



6 V 4AT_ 

In addition, when S(2p(j)) = 0, we have from Eqs. (37-42) 

iA - = ln2+ 6^ + 60^ + °(^)- <™ 

This indicates that the magnitude of the positive MC is gradually increased (e.g., the saturation value of AG is raised 
above 1) when contributions from higher-order terms (i.e., terms oc 1/N with k > 1) are included, though they are 
negligibly small. 



D. Discussion 

Our results for the MC are in good agreement with experimental measuremeuts.lliS For instance, a positive MC, 
is observed in the VRH regime of both macroscopically large I^Oa-^ samplesEEl and compensated n-type CdSe.E2l 
Moreover, saturation in the MC as the field is increased is also reported in Ref. 20. 

The results for AG(^ presented in this work are consistent with theoretical studies based on an independent- 
directed-path formalisma and a random matrix theory of the transition strengths .u The advantages of our results 
include: (i) they provide explicit expressions for the first two dominant contribution to the MC, as a function of the 
magnetic field; (ii) they provide straightforward determination of the period of the oscillation of the MC; (iii) they 
provide explicit scaling behaviors (i.e., the dependence on the hopping length and the orientation and strength of the 
field) of the low- flux MC in quasi 1-D, 2D and 3D systems; and (iv) they allow us to make quantitative comparison 
with experimental data. Finally, it is important to emphasize that our analytic results [Eqs. (48), (49) and (52)] for 
the MC are equally applicable to any dimension, since the essential ingredient in our expressions is the QI quantity 
S^ r \ which takes into account the dimensionality. 

In appendix B, we outline the computational scheme using the second model of disorder, i.e., e; is randomly chosen 
from the interval [-W/2, W/2]. The moments obtained in this case are the same as those presented in Eqs. (37) and 
(38). Therefore, the result for the MC remains unchanged. 
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IV. QUANTUM INTERFERENCE AND THE SMALL-FIELD MAGNETOCONDUCTANCE ON A 

THREE-DIMENSIONAL CUBIC LATTICE 



A. Sums over forward-scattering paths 



Let 

S m ,n,i (— in 3D) be the sum over all phase factors associated with directed paths of m + n + Z (— r) steps 
along which an electron may hop from (0,0,0) to the site (to, n, I). Again we assume to, n, and I > 0. In other words, 
electrons can now also hop in the positive z direction. The vector potential of a general magnetic field (B x , B y , B z ) 
can be written as 

A = ^(zB y - yB z , xB z - zB x , yB x - xB y ). 

Also, a/2n, b/2ir and c/2ir represent the three fluxes through the respective elementary plaquettes on the yz-, zx- 
and xy- planes. To compute 5 TO ,n,;, we start from the following recursion relation 

S m ,nJ = ^ ^ ^p,g,^TO,rt,i exp ii J S p>q> l-1, (55) 

where j4p, g ,i— >m,n,i is the sum over all directed paths starting from (p,q,l) and ending at (m,n,l). The physical 
meaning of Eq. (55) is as follows. The site (m, n, I) is reached by taking one step from (p, q,l — l) to (p, q, I), acquiring 
the phase i(qa — pb)/2, then traversing from (p,q,l) to (m,n,l) on the z = I plane. After some calculation, we find 
that 



^- P , q ,l-^m,n,l CXp \ I 



(to — p)(lb — qc) + (n — q)(pc — la) 



^m— P ,n—q{p)i 



(56) 



where S m - P ,n- q (c) is defined as shown in Eq. (10). By applying Eq. (56) I times, we obtain a general formula of 
Sm,n,i for m, n, I > 1 in terms of the fluxes a, b and c as 



S m ,n,i(a,b,c) = exp 



nla + Irab + mnc 



C m ,n,i(a, b, c), 



(57) 



where 

£ra,n,l(«>M) 



Pj+i Qj+i 



L pi,qi(c), (58) 



H ex P ^fe a + ( m - Pj) b + Pj(lj+i - Qj)c]} L Pj+1 - Pj , qj+1 - qj (c) 
Pj =o qj =o 

with pi + i = to, qi + i = n, and the L Pj9 (c)'s are defined as in Eq. (12). 

It is clear that S m , n ,o — S m , n (c), S m .o,i = S mi i(b), and So, n ,i = S n ,i{a). Also, the following symmetries hold: 

S m ,n,i( a > & > c ) = a > c ) = 5 m ,i,n(a, c, 6) = <S;, m , n (&, c, a) = Sj,„, m (c, 6, a) = S„,j, m (c, a, 6). (59) 

When there is no magnetic flux, 

^(0,0,0)= ^ + "+/^ ^ 
to! n\ l\ 

gives the total number of (to + n + Z)-step paths connecting (0, 0, 0) and (to, n, I). 

We have obtained explicit expressions for many 5 mi „ ; ;, and here we explicitly present only the first few S, since S 
have long expressions for larger to, n, and I. 



<5i,i,i 

$2,1,1 



= 2 



a+b—c b+c—a c+a—b 
cos h cos h cos 



2 [cos — + cos ( — — b 
12 \2 



cos ( c) + cos ( — ± b =F c 



52,2,1 = 2 + 2 cos a + 4 cos(a - 6) + 2 cos(a - c) + 2 cos(6 - c) + 2 cos(a - 2c) + 2 cos(6 - 2c) 

+ 2cos(a -b±c) + 2cos(a + b - 2c) + 2cos(a -b±2c), 

6 + [4 cos a + 2 cos 2a] + J^ (2) [4 cos(a - (i) + 4 cos 2 (a - (3) + 2 cos(a - 2(3) + 2 cos(2a - (3)] 

+ 2^ (3) [cos(a + /3- 7 ) + cos 2(a + (3 — 7) + cos(a + /3 — 27) + cos(a ± 2/5 =F 27)] . 



'2,2,2 



Here denote sums over a = a, 6, c; (a /3) = (ab), (fee), (ca); and (a/37) — (afec), (oca), (cab); for i = 1, 2, and 3 
respectively; and, for instance, the term cos(a/2 ± b =F c) means cos(a/2 + 6 — c) + cos(a/2 — 6 + c). 
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B. Low- flux limit 

In the very-low-flux limit, and calculated exactly to second-order in the flux, we obtain the logarithm of 5 m>ri> ;, the 
3D analog of the harmonic shrinkage of the wave function, as 

\nS m ,n.i = In J\f — — [nla 2 + Imb 2 + mnc 2 + m(lb — nc) 2 + n(mc — la) 2 + l(na — mb) 2 ~\ . (60) 

This generalizes the 2D harmonic-shrinkage of the wave function obtained in Eq. (13). When m = n = I = r/3, we 
have 

ln<S ro , ro>m = In - |r 2 (a 2 + b 2 + c 2 ) + j [(b - c) 2 + (c - a) 2 + (a - b) 2 } \ . (61) 

These results generalize to 3D the 2D results obtained in Sec. II B. 

C. Interference patterns on diagonal sites 

In order to see how the interference patterns vary according to the orientation of the applied field, we focus on 
S m .m,m (i.e., S on the body diagonals). We now examine two special cases: B|| = B || (1,1,1) = (</>, (j>, 4>)/2tt 
and = B _L (1,1,1) = (<fi/2, (f>/2, — <fi)/2ir, namely fields parallel and perpendicular to the (1,1,1) direction, 

respectively. Their 5 mjlriim 's are denoted respectively by l\ m and T^ m and have been computed to high orders. Here 
we only present the first few: 

jj = 6 cos 0, 

j| = 36 + 42 cos 20 + 12 cos 40, 

4' = 864 cos + 528 cos 36* + 216 cos 50 + 54 cos 70 + 18 cos 90, 

j| 2 = 7308+ 12504 cos 20 + 8082 cos 40 + 4032 cos 60 + 1740 cos 80 
+ 672 cos 100 + 216 cos 120 + 72 cos 140 + 24 cos 160; 

and 

l£- =4cos0 + 2cos20, 

l£- = 14 + 12 cos + 16 cos 20 + 12 cos 30 + 12 cos 40 + 8 cos 50 + 10 cos 60 + 4 cos 70 + 2 cos 80, 
Ig = 76 + 204 cos + 176 cos 20+ 180 cos 30+ 156 cos 40 + 156 cos 50 + 136 cos 60 
+ 128 cos 70 + 102 cos 80 + 84 cos 90 + 68 cos 100 + 64 cos 110 + 48 cos 120 
+ 40 cos 130 + 26 cos 140 + 20 cos 150+10 cos 160 + 4 cos 170 + 2 cos 180, 
1^2 = 1372 + 2464 cos + 2606 cos 20 + 2420 cos 30 + 2502 cos 40 + 2288 cos 50 + 2288 cos 60 + 2068 cos 70 
+ 2046 cos 80 + 1788 cos 90 + 1758 cos 100 + 1532 cos 110 + 1498 cos 120 + 1264 cos 130 
+ 1174 cos 140 + 964 cos 150 + 894 cos 160 + 724 cos 170 + 642 cos 180 + 512 cos 190 
+ 450 cos 200 + 340 cos 210 + 296 cos 220 + 228 cos 230 + 178 cos 240 + 128 cos 250 
+ 94 cos 260 + 56 cos 270 + 40 cos 280 + 20 cos 290+ 10 cos 300 + 4 cos 310 + 2 cos 320. 

where = ^/2. It can be seen that xj m and T^ m exhibit quite different behaviors as shown in Fig. 7 where we plot l\ 

through l\ 2 and 1% through 2^. Notice that the period in <j> for zj m is 2ir (4ir) for even (odd) to, while the period 
for T^ m is Att for any to. Therefore, the periodicity for the MC in 3D is identical to that in 2D. 

We have also computed x| m and X^ m (to = 1,2,..., 300) for <f)/2it = 3/5 and (\/5 — l)/2 and find that their 
behaviors are insensitive to the commensurability of 0, unlike the case on a square lattice. Physically, this can be 
understood because two randomly chosen paths have a higher probability of crossing (and thus interfering) in 2D 
than in 3D; thus making QI effects less proBajiHced in 3D than in 2D. A similar situation occurs classically (e.g., 
multiply-scattered light in a random mediumMi 16 !) . 
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D. Small-field magnetoconductance 



For a 3D system, the relative MC, AG(a, b, c), now reads 

AG(a,b,c) 



S mtTl j(2a, 2b, 2c) 



1 - 



J\f 



(62) 



The above general expression is valid for any ending site as well as arbitrary orientation and strength of the magnetic 
field. From Eq. (61), in the small-field limit and at ending site (r/3, r/3, r/3), we have 



AG 



1 



_ v 'r 2 (a 2 + b 2 + c 2 ) + — [(b - c) 2 + (c - a) 2 + (a - &) 2 ] 



(63) 



which is applicable for any orientation of the field. Below we focus on two special orientations of the field: and 
B„. 

For very small </>, we have from Eq. (61) 

hiZ^ ~ In 

and 



[(r/3)!] 3 144 



r 2 (r + 1) , 



lnTjJ ~ In ■ 



! 1 

[(r/*W ~ 72 

The 3D behavior of the low-flux MC thus becomes clear: for B i 



2 J.2 



AG{<j>) 



^r 3 / 2 0, 



and for B 



AG(0)~~r«£. 

These results can be interpreted as follows: the effective area exposed to B^ is larger, 

Af ~ r 3 / 2 , 

similar to the 2D case where AG ((f)) oc r 3 / 2 <ft; while the effective area Ajj ff exposed to B|| is smaller, 



(64) 



(65) 



(66) 



(67) 



Af - r, 

thus closer to our quasi-lD ladder case with AG(4>) ocrcf). 

As a numerical test of Eqs. (66) and (67), in Fig. 8 we show AG calculated directly from Eq. (62), versus r 3 / 2 <fi 
in (a) and versus rtfi in (b), respectively for several values of B^ and B|| with hopping length r = 3, 6, . . . , 600. The 
collapse of all the data into a straight line verifies the scaling of the low-flux AG presented above. 



E. Average of the magnetoconductance over angles 

In a macroscopic sample, the conductance may be determined by a few (instead of one, as considered before) critical 
hopping events. As a result of this, the observed MC of the whole sample should be the average of the MC associated 
with these critical hops. Thus, in 3D systems it is also important to take into account the randomness of the angles 
between the hopping direction and the orientation of the applied magnetic field.Ej 

To theoretically investigate the effect of the average over angles on the MC, we consider all possible relative hopping 
directions with respect to that of the magnetic field, or equivalently, the continuously-varing orientation of the field 
with respect to a fixed hopping direction. We adopt the latter below: the ending site of all hopping events (with 
the same hopping length r) is located at the diagonal point (r/3, r/3, r/3) and the magnetic field can be adjusted 
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between the parallel and perpendicular directions with respect to the vector d = (1, 1, 1). Our interest here is in the 
MC averaged over angles, denoted by AG, in the low-field limit. Recall that the magnetic field is B = (a, b, c)/2tt, 
and from Eq. (63), we have 

AG = r B \fl + r sin 2 lu , (68) 
3\/3 



where B = \J a 2 + b 2 + c 2 /2tt is the magnitude of the field and lu is the angle between B and d. By averaging over 
the angle to, we obtain 

7T/2 



4 r*l* , 

AG(B) = —= rB I Vl + r sin 2 to duo 
3v3 Jo 



4 



3^3 \2' y/F+TJ ' 



Vf+1BE , (69) 



where E(n/2, \/r/\/r + 1) is the complete elliptic integral of the second kind. When r is large, E(w/2, \frj\fr + 1) ~ 1 
and we therefore have 

AG(B)^r 3 / 2 B. (70) 

Equation (70) means that the dominant contribution to the MC stems from the critical hop which is perpendicular 
to the field. This is understandable through our earlier observation that the effective area enclosed by the electron 
is largest when B is perpendicular to d. From the above analysis, we conclude that in 3D macroscopic samples the 
low-field MC should in principle behave as r 3 / 2 B. 



V. CONCLUDING REMARKS AND SUMMARY OF RESULTS 



In closing we briefly address four issues. First, although relevant measurable quantities such as ISVn.nl 2 and |iS m ,n,z| 2 
are gauge-invariant, the transition amplitudes are gauge dependent. As an illustration, the transition amplitude will 
be L m , n [Eq. (12)] if we use the Landau gauge A = (0, Bx) on a square lattice. The notation S m , n {L m ,n) refers the 
use of the Symmetric (Landau) gauge. Similarly, the transition amplitude will read C m ^ n j [Eq. (58)] if we use the 
gauge A = (B y z, B z x, B x y) on a cubic lattice. 

Second, returns t^ithe origin (see, e.g., Refs. 1 and 22-24) become important for less strongly localized electrons, 
and their QI effectsESO can be incorporated in our approach. 

Third, the main limitations of our study in the case with impurities are the following: no inclusion of spin-orbit 
scattering effects (for this see, e.g., Refs. 7-9 and references therein), and no explicit inclusion of the correlations 
between crossing paths, as discussed in Refs. 4 and 7. However, these correlations are negligible when spin-orbit 
scattering is present .13 

Fourth, besides analytical closed-form results in 2D, this work presents exact results for 3D systems, e.g., AG = 
(27r/3\/3)rS(l + rsin 2 to) 1 / 2 [Eq. (68)]. These results can provide further tests of the quantum interference effects. 
This can be done by measuring the MC of bulk samples (which are small enough that only a single critical hop 
is allowed) in various orientations of the field. By doing this, one can then determine the values of r and to (and, 
hence, also the direction of this critical hop). Therefore, the small- field behaviors of the MC with fields parallel and 
perpendicular to the direction of this critical hop can be measured and compared to our predictions. This could 
potentially be very useful. 

In summary, we present an investigation of quantum interference phenomena and the magnetic-field effects on the 
MC resulting from sums over directed paths on resistor networks in 2D and 3D. The principal results include: (1) 
an exact and explicit closed-form expression for the sum over forward-scattering paths S^ r ' to any site on a square 
lattice, which is the essential QI quantity in both uniform and disordered cases, (2) an explicit formula for for 
electrons hopping on a cubic lattice, (3) the low- flux behaviors of in both 2D and 3D, (4) the exact summation 
of the dominant winding paths in 2D, (5) compact, analytic results for the positive MC as explicit functions of the 
magnetic flux which are valid in any dimension, (6) the small-field behaviors of the MC in quasi-lD, 2D and 3D, and 
(7) an analytic result for the small-field MC in 3D incorporating the randomness in the relative angles between the 
hops and the applied field. They provide analytical and explicit closed-form results concerning the hopping transport 
of strongly localized electrons subject to an external magnetic field in the macroscopic regime. We hope that our 
results stimulate further work (e.g., inclusion of spin-orbit effects on lattice path integrals) on exact results in 2D and 
3D systems. 
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APPENDIX A: DERIVATION OF THE MOMENTS (Jf? iN {<f>)) AND ( | J m ,n(4>)\' 2P ) 

In this appendix, we outline the derivation of the moments shown in Eqs. (32), (33), (35), and (36). First, let us 
derive Eqs. (32) and (33). Let Sj = exp(i<&j), where <&i is the sum over all phases along path i and i = 1,2,..., AT. 
Note that for every path, the probability for obtaining the overall phase factor ±s t is P±. Now 

J 2 = 

and 

|J|2= y- / I----/- 

\l=l / \l=l / 

where ji = ±1 with probability P±. If the number of ji = —1 is k, the overall probability is P+~ P 1 ^ and there 
are combinations among ji (i = 1, • • • , N). For k = and k = N, there is only one combination producing 
S(2<j>) + 2^2 i: £j Si Sj for J 2 and also only one combination producing N + Si/sj for |J| 2 . When 1 < k < 

N - 1, for J 2 there are V_ = 2C"j 2 combinations producing S{2(f) - and N + = Cj? - 2C^ 2 

combinations producing S(2<fi) + 2J2i^j s iSj. Also, when 1 < k < N — 1, for |J| 2 there are 7V_ combinations 
producing N — s i I s j an d N + combinations producing N + s i I s j ■ In N- , the factor 2 comes from the two 
possible ways: +Sj — Sj and — Sj + Sj. C^~ 2 comes from arranging (k — l)'s minus signs among the (N — 2) s;'s left 
(/ + i + j). Therefore, the overall average is V = Pf + P N + J2k=i( N + - N -) P +~ k P- = (p - v) 2r . We thus have 

( J 2 ) = S(2cj}) + V \ 2y St s 




where we have used the relation J2k=o C k p +~ k p - = ( p + + P -) N = L Similarly 



By exploiting 



■ ■ S 3 



2j2s iSj =S 2 (ct>)-S(2ci>) 



and 

we thus obtain Eqs. (32) and (33). 

For /j, = v = 1/2, (namely, each site contributes either +1 or —1 with equal probability), P± = 1/2 for every path 
i. The total number of sites that can be visited is T — (m+ l)(n + 1) — 1 (the initial site is not counted). Therefore, 
the total number of all possible impurity configurations is 2 T . Let us now focus on the 7$ Sj with i = 1, 2, • • • , N. The 
total number of configuration sets is 2 N since each 7, can be either +1 or —1 with equal probability. Among the 2 T 
impurity configurations, there are always 2 T /2 N producing a set of 7^ for every possible set of 7^. For instance, 
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T = 3 and N = 2 in the simplest case S\^. Among 2 3 = 8 impurity configurations, two of them give +s\ + s 2 , 
+si — s 2 , —si + S2, and — si — s 2 . We thus have 



P / N 



{ 7i } L \*=1 / \i=l 



and similarly, 



{ 7 ,} L \i=l / \i=l 1 



The summation in {7^} is over all possible configuration sets of 73. In other words, the average of A over disorder 
(i.e., (A)) means to sum over all possible ji = ±1 for the desired quantity A, and then divide the sum by 2 N . We 
thus obtain 



u 4 > = w) + t|^£^ 

(j8> = 5(80) + 6T2jE s ^' + ^E s ^l + ir|j)2 E a * a i a * + w. ^ a< w/ - 

U 10 > = W) + ^E-?«? + ^E« + eiW.E W + (Si .£ 



E 



2222 

sf si si sf 



+ 



10! 



E 4 2 2 2 
s t s j s k s l 



10! 

(2P 



E 



22222 

S i S j S fc S Z S m i 



_ 4 4 2 

J afe (41)2 21 ^ Sl S i Sk 



and 



<|J|*>=tf(2tf-l) + £* 

, s 2 

( I J| 6 ) = N (6N 2 - 9N + 4) + 3 (37V - 4) ^ -\ , 

2 

( I J| 8 ) = AT (24iV 3 - 72iV 2 + 82N - 33) + 4 (18iV 2 - 57iV + 49) ^ 4 + E ' 



+ 6E 



<j 4 s 2 s? 



E s fc s ; 
s 2 s 2 ' 



( I J| 10 ) = N (120iV 4 - 6007V 3 + 12507V 2 - 1225iV + 456) + 20 (30iV 3 - 1657V 2 + 3257V - 224) ^ % 



+ 5 (57V - 8) % + 10 ( 15N - 32 ) E 



„2 2 
S j s k 



S 4 



2 2 4 



300 (3JV - 8) ^ 



s • 

s 2 s 2 
s 2 s 2 

i^j^k=£l k 1 



Noticing that J^iLi S T — J2iLi ^l s T = S(m<f)), it can be derived that: 



£ S ^ 2 = V^)-^)], 



E s * s2 = 5(2^)5(4^) -5(60), 
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E^4 = ^ 2 (40)-5(80)], 
E^ 2 = S(20) 5(60) -5(80), 
E*!* 2 = 5(20) 5(80) -5(100), 
E*M = 5(40) 5(60) -5(100), 

E * 2 a * * 2 = ^ s 3 ( 2< « - ^ 5 ( 2 ^) w) + 1 w). 

E s * s ) s l = \ s2 ^) - - \ s2 (^) + 5 ( 8 ^)' 

E s i s2 s l = \ s2 ^) - s w) - \ + 5 ( 10 ^)' 

i^j^k 



E s * 8 i s2 = \ S W) &W>) - \ 5(20) 5(80) - 5(40) 5(60) + 5(100), 

i^j^k 

E s2 s2 s I s2 = ta s4 ^ i s2 ^ s ^ + \ s <w + 1 s2 ^ ~ 1 w)> 



E 4 s 2 s l^ = \ S 3 (20) 5(40) - \ 5 2 (20) 5(60) - \ 5(20) 5 2 (40), 



+5(20) 5(80) + | 5(40) 5(60) - 5(100), 

E s2 s l * 2 * 2 s ™ - 55(2< ^ ~ h. 53(2<w sm + 1 s2m S{ ^ ] + 1 sm s2m 

- \ 5(20) 5(80) - \ 5(40) 5(60) + \ 5(100), 

5? 



E^- 5 2 (20)-iV, 

e4=^ 2 ( 4 <a)-^ 



E 4V= E #-|5 2 (20)5(40)-5 2 (20)-i5 2 (40) + 7V, 

2 2 

E S^i = I 54 ( 2 ^ - ( N ~ 2 ) - J S 2 (20) W) + I S " + ^ - 3). 

By utilizing the above relations, we have obtained the results for the moments presented in Eqs. (35.1)-(35.4) and 
(36.1)-(36.4). 



APPENDIX B: A DIFFERENT MODEL OF DISORDER 



In this appendix, we study another model of diagonal disorder: ej uniformly distributed between — W/2 and W/2. 
Our focus is on the analytical computation of the leading terms (terms oc N p ) in the moments (|J(0)| 2p ) and 
( | J(0)| 2p ). Indeed, the scheme presented below applies equally to our first model of disorder, where ej can take two 
values, +W and — W, with equal probability. 



20 



Let us write 



J (<W = E^ Si ' 



where for each path i 



*-in-£V 



i=i 



It is reasonable to choose rji from a Gaussian distribution of zero mean and unit standard deviation. We therefore 
have {rjj n ) = 1 and {rft n+1 ) = 0. For <f> = 0, we have 



(AT \ 2 P 2p / N \ 
E* =n EM 
*=1 / fe=l \i*=l / 



It is found that the leading term in ( | J(0)| 2p ) comes from all terms having the form Yii=i Vi an( i there are Cp distinct 
terms of this type, each term has a coefficient (2p)\/2 p . Therefore, we obtain for the leading term 



<k(0)| 2p H%^ = (2 P -l)!!^. 



For <j) ^ 0, we have now 



ttIYv Uf r 



fe=l 



The contributions to the leading terms involve different factors, as shown below. There are Cp terms like 

(eH OH)' 



and they contribute 



There are C^_ 2 2 terms like 



MP 

{p\f ^=p\N*>. 



and they contribute 



/ p-i \ 

\ i=2 / 



iP tt 

S l' i=2 V 



2 J (p-2)! (l!) 2 ^ P '(2.1!) 27V 6 (20) ' 



Similarly, there are C^ 4 4 terms like 



\ / x i f- 2 i \ 
vl 4 r& s l II ^ Si ) | ^i' 4 II Vi 7, J ' 



i=3 



2' J=3 " J 
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and they contribute 



p\\ 2 NP- 4 S 4 {2cj)) 4!Cf 



2 2 y (p-4)! (2!) 2 



(2 2 -2!) ; 



•iV p - 4 S 4 (2c/>). 



In general, there are C^_ 2 ^, k terms of the following type 



■ y p-2k 



p—k 



n in? n ^ > 

\i=l i=fe+l / \i=l *' i=fc+l V 
each one with a coefficient (p!/2' c ) 2 . The contribution to the moment is therefore 



pi 



2 k J {p-2k)\ {k\f 
Notice that we have utilized the fact that 



(2 fc fc!) 2 



k s 2 



All difFcront and l t \i=l '« / 

always contains S 2k {2(f)) / (k\) 2 . Totally, we thus obtain for the leading terms 



vfe = 



(2fc)!C 2 P ; 
(2 fe fc!) 2 



2/,- 



S{2<j>) 
N 



2k" 
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FIGURE CAPTIONS 



FIG. 1. (a) Starting from (0,0) on a square lattice, for forward-scattering paths of four steps, electrons can end at five 
sites: (4,0), (3, 1), (2,2), (1,3), and (0,4). Their corresponding S m ,n are also shown. The arrows specify the electron hopping 
directions (only moving to the right and upward are allowed in the directed-path model). Notice that the symmetry S m ,n = S n ,m 
holds. 52,2 has the strongest interference among them because the number of paths ending at (2, 2), and the area they enclose, 
are both the largest, (b) Six different directed paths connecting (0,0) and (2,2) and their separate phase-factor contributions 
to S 2 ,2; the total equals 1 + 1 + e l4> + e _i * + e 2 ^ + e~ 2 ^ = 2 + 2 cos + 2 cos 20. 



FIG. 2. Plot of m versus 0/2-7T (denoted by short bars in order to visualize them better), between and 1, such that 
l2m{4>) ~ 0; for m — 1,2, ... , 20. Note that the smallest one is always l/2m and the number of zeros increases rapidly when m 
becomes larger. The properties of lira described in Eqs. (19-22) are exhibited in the figure. For instance, when <j)/2n — 1/5, 
hm — for m = 3 + 5n and m = 4 + 5n with n > (namely, m = 3,4, 8, 9, 13, 14, 18, 19, . . .). 



FIG. 3. hm for various m as functions of the flux through each elementary plaquette, $ = 4>/2n in their respective full 
period. Notice the 2n (4tt) periodicity in <f> for even (odd) m. In (a), we plot I2, 14, ■ ■ ■ , I12, lis, and ho- To show the behavior 
of the rapid small- magnitude fluctuations around zero of hm(&) for ^7 < $ < | when m is even and for ^7 < $ < when 
m is odd: In (b), we plot 7io (top), lis (middle), and /3s (bottom) for $ in their respective interval [^, §]. In (c), we plot 7i2 
(top), I20 (middle), and I40 (bottom) for $ in their respective interval [^-, Only some restricted ranges in the vertical 

axes are exhibited. From these figures, we clearly see the general properties for the behavior of 7 2m described in Sec. II C. 



FIG. 4. P m ,m (for m = 1, 2, . . . , 6) as functions of the flux through each elementary plaquette, $ = (j>/2n. 

FIG. 5. The relative magnetoconductance AG(<j>) versus <j>/2n for hopping between (0,0) and (r/2,r/2) for several system 
sizes. From (a) to (d), the hopping length r corresponds to 4, 10, 20, and 50, respectively. Inserts show AG((j>) for <j> between 
and the corresponding saturated field 4>/2ir = l/2r. It is observed that for large systems (i.e., r large), AG((j>) rapidly 
approaches the saturation value 1 even at 4>j2-K, which is less than l/2r. 

FIG. 6. (a) AG versus r 3/2 B in 2D with the hopping length r = 2, 4, ... , 1000, and (b) AG versus rB in quasi-lD with 
r — 2, 3, ... , 500, for various small values of B. All the data nicely collapse into a straight line, which verifies the scaling 
behavior of the small-field AG: (\/3/6)r 3/2 in 2D and ( v / 3/3)r0 for quasi- ID systems. The distance between these data and 
the solid reference line reflects the prefactor: ^3/6 in (a) and a/3/3 in (b). 

FIG. 7. Sums over forward-scattering paths between two diagonally-separated sites on a 3D cubic lattice: ij through l\ 2 
for B = (<f>,(f>,(f>) and 1$ through f° r B = (4>/2, 4>/2, ~(f>), as functions of 4>/2tv. Note that while lj m has the 27r (4-7t) 
periodicity for even (odd) m, 1^ always has a period 4-7T. 

FIG. 8. (a) AG versus r 3/2 B for Bi = B (1/2, 1/2, -1), and (b) AG versus rB for By = B (1, 1, 1) for several values of 
B and hopping length r = 3, 6, ... , 600. The collapse of all the data into a straight line verifies the scaling of small-field 
AG : (\/2/6)r 3//2 for Bi and (l/3)r0 for By. The distance between these data and the solid reference line reflects the 
prefactor: \72/6 in (a) and 1/3 in (b). 
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